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Abstract 

One-dimensional Yang-Mills Equations are considered from a point of view of a class 
■ of nonlinear Klein-Gordon-Fock models. The case of self-dual Nahm equations and non- 

self-dual models are discussed. A quasiclassical quantization of the models is performed 
by means of generalized zeta-function and its representation in terms of a Green function 
diagonal for a heat equation with the correspondent potential. It is used to evaluate the 
. functional integral and quantum corrections to mass in the quasiclassical approximation. 

Quantum corrections to a few periodic (and kink) solutions of the Nahm as a particular 
case of the Ginzburg-Landau (phi-in-quadro) and Sin-Gordon models are evaluated in 
arbitrary dimensions. The Green function diagonal for heat equation with a finite-gap 
potential is constructed by universal description via solutions of Hermit equation. An 
alternative approach based on Baker- Akhiezer functions for KP equation is proposed. The 
generalized zeta-function is studied in both forms; its derivative at zero point, expressed 



in terms of elliptic integrals is proportional to the quantum corrections to mass. 

> : 

1 Introduction. 

1.1 General remarks. 

We consider one-dimensional field theory, based on nonlinear Klein-Gordon equations, arising, 
for example of Sine-Gordon (SG) case, in kink models for crystal structure dislocations [1] or in 
a context of a relativistic models [2]. Popular operator constructions of quantum field theory 
such as Yang-Mills one have reductions to one dimensional models [3]. Some kind of embedding 
of such theory into the multidimensional one is possible: Atiyah-Drinfeld-Hitchin-Manin-Nahm 
construction may appear as an equivalence between two sets of self-dual equations, one as 
described above in one dimension, the other in three dimensions (reduced from a Euclidean 
four dimensional theory by deleting dependence on a single variable) [I]. 

A class of nonlinear Klein-Gordon equations in the case of static one-dimensional solutions 
is reduced to 

(j," -V'{(j>) =0,<p = <p(x),x e R. (1.1) 
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Suppose the potential V(4>) is twice continuously differentiable; it guarantees existence and 
uniqueness of the equations correspondent to (II. ip Cauchy problem solution. The first integral 
of the equation (11.11) is given by 

W = ^' 2 -V^), (1.2) 

where W is the integration constant. The equation (ll.2p is ordinary first-order differential 
equation with separated variables. As the phase method shows, solutions of this equations 
belong to the following families: constant, periodic, separatrix and the so-called "passing" one 

nu. 

In the case of Nahm equation 

^v(0) = y; w = ~y> (L3) 

and for Ginzburg-Landau the "potential" is 

W) = f(0 2 -y) 2 , (1.4) 
while in the case of Sin-Gordon (SG) model it is 

VsM = ^ + cos(-\lfx)). (1.5) 

We modified this last model to fit it with the first one at small values of the constant g, namely 

13m 4 . 0x 

Vs 9 = ^ + V gl + 0(g 2 ). 

1.2 One-dimensional reduction of Yang-Mills theory. 

Starting with the full Yang-Mills equations in four Euclidean dimensions, 

DpT^ = 0, (1.6) 

for the gauge fields T M = T+, where 

= T v<ll - Tp, v - i[T M , T„], = d„- i[T M , $]. 

and demanding all fields be independent of three of the variables xj., k=l,2,3 setting x^ = z, 

&=[T j [T j ,T k ]\, [T k ,%] = 0. (1.7) 
One cam easily check that the self-dual equations, 

-j^ = ±£ijkTjT k , (1.8) 
corresponding to the model, imply Eqs. (II. 7p . Starting from the simplest solution of the system 

(USD 

Tt = Mi (1-9) 



2 



built on a base of Pauli matrices one arrives at the Euler system for (j>i{y) that is solved in 
Jacobi functions. The solutions are dressed by the gauge-Darboux transformations [BJ. 
For the second order equations the ansatz similar to (II. 9p 

Ti — (f)(z)a>i,i — 1,2,3, (1.10) 

with a constant matrices ctj and a convenient choice of scaling leads to the pair of equations 




(1.11) 



The second order equation for <fi(z) enters as a particular case into the class of nonlinear (11.11) 
with V'(<f)) = 2(p 3 , or into Lagrangian of GL model with V from (11.31) . 



1.3 Feynmann quantization of a classical field. 

An attention to Feynmann quantization formalism of a classical field was recently attracted in 
connection with a link to a SUSY quasiclassic quantization condition [8] suggested in [9]; see, 
however, [10] . 

Historically, quantum corrections started from a [12], see also [13] ■ Important development 
concerns the Jacobi variety structure [13J. 

In the papers of V.Konoplich [15] quantum corrections to a few classical solutions by means 
of Riemann zeta-function are calculated in dimensions d > 1. Most interesting of them are the 
corrections to the kink - the separatrix solution of the field 4 (GL) model [IB] . 

The method of [15J is rather complicated and it is desired to simplify it, that was the main 
target of our previous note [IT]. We applied the Darboux transformations technique with some 
nontrivial details missed in |15j . 

The suggested approach open new possibilities; for example it allows to calculate the quan- 
tum corrections to matrix models of similar structure , Q-balls [H] and periodic solutions of 
the models. The last problem is posed in the review |llj . 

The approximate quantum corrections to the solutions of the equation (11.11) are obtained 
via the Feynmann functional integral method evaluated by the stationary phase analog [19]. It 
gives the following relation 

^- S ir ] -7m- < 112 > 

where S qu denotes quantum action, corresponding the potential V(4>), A - some quantity de- 
termined by the vacuum state at V{<p) = 0, and detD is the determinant of the operator 

D = -%-Ay + V"{<Kx)). (1.13) 

The argument y G R d ~ x stands for the transverse variables on which the solution 0(x) does not 
depend. The operator D appears while the second variational derivative of the quantum action 
functional (which enter the Feynmann trajectory integral) is evaluated. For the vacuum action 
S vac the relation of the form (11.121) is valid if S qu is changed to S vac and D is replaced by the 
"vacuum state" operator D = —d 2 x — A y + v. Then, the quantum correction 

AS qu = Sq U S vac , (1.14) 
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is obtained by the mentioned twice use of the formula (I1.12p as 



H , . det D 

A5 '"=2 ln W' (115) 

Hence, the problem of determination of the quantum correction is reduced to one of evaluation 
of the determinants D and Dq ratio. The methodic of the evaluation will be presented in the 
following section. 



1.4 The generalized Riemann zeta-function and Green function of 
heat equation. 

The generalized zeta-function appears in many problems of quantum mechanics and quantum 
field theories which use the Lagrangian £ = (<90) 2 /2 — V((f>) and it is necessary to calculate a 
Feynmann functional integral in the quasiclassical approximation. 

The scheme is following. Let the set {A„} = S be a spectrum of a linear operator L, 

ln(detL) = lnA ™' ( L16 ) 

\ n es 

where the sum in the r.h.s. is formal one. 

The generalized Riemann zeta-function (l(s) is defined by the equality 



Cl(s) = V- (1.17) 



A„eS 

This definition should be interpreted as analytic continuation to the complex plane of s 
from the half plane Res > a > in which the sum in (11.171) converges. Differentiating the 
relation (11 . 17)) with respect to s at the point s = yields 

ln(detL) = — Ci(0). (1.18) 

The generalized zeta-function (13.891) admits the representation via the diagonal Ql of a Green 
function of the operator d t + L. 

A link to the diagonal Green function (heat kernel formalism) has been used in quantum 
theory since works by Fock ([21]). There is a representation in terms of the formal sum over 
the spectrum 

g L (t,r,T ) = ^exp[-A n t]^ n (rX(r ), (1.19) 

n 

where the normalized eigenfunctions ip n (r) correspond to eigenvalues A n of the operator L. 
Let us next define 

1L (t) = E fe e~ Xkt = $9L {t, r, r) dv (1.20) 
The Mellin transformation yields the generalized zeta function of the operator L: 

i r+oo 

Us) = m f- l lL (t)dt. (1.21) 
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Returning to the operators fll~T3]) L -> D,D : D = -A + ^"(0 o (x)) and D = - = -£5 - 
A y + v, we pose the main problem for a periodic potential u(x): 

— + D l \g D (t, x, x ) = 6(t)5(x- x ) (1.22) 

where D\ = —£2+ u(x), V < V m , that means a necessary presence of a continuous part of 
the spectrum A G [Ao, +00). The basic relation (11.151) points to a necessity of evaluation of the 
determinants of the operators 

D = D + u(x), D = -dl-A y + v (1.23) 

where A is a positive number and x G R is one of variables, while y G Rd-i is a set of other 
variables. The operator A y is the Laplace operator in d-1 dimensions, u(x) is one-dimensional 
potential that is defined by the condition 

V"(Mx)) = ^ + u(x), (1.24) 

where <fto(x) is the classical static solution of the equation of motion. 
The Green function for a Hermitian operator Df = Di is : 

g Dl (t, x, x ) = ^ e ~ Xkt i>k (x) i>l (^0) © (t) ■ 

k 

The generalized zeta- function, defined by the relations (ll.20|ll.2ip . will be referred as the 
zeta-function of the operator D. 

From the relation (11.201) for the function 7d(£) it follows an important property of multi- 
plicity: if the operator D is a sum of two differential operators D = D\ + D2, which depend on 
different variables, the following equality holds 

J D {t) = j Dl (thD 2 (t)- (1-25) 
Generally, for a constant potential v (see Sec. 3), for the one- dimensional case, 



We will need the value of the function 7£>(i) for the vacuum state, when the operator D = —A y 
is equal to the d-l-dimensional Laplacian. In this case v — 0, hence 

lD (t) = t^tt / rfkexp(-|k| 2 t) = {Avt)-^l\ (1.27) 

Combining as in ( 11 .251) yields 

InM = ^# (1.28) 

2v Tit 



Then, 

1 p+00 Y(s -d) 1 

CdM = W)1 f ^ a ' m = ■ (L29) 

A quantum correction to the action in one-loop approximation for the classical solution 4>(x) 
is calculated via zeta-function by the formula 

A5 fltt = -?[(l,(0)-(^(0)]/2. (1.30) 
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2 The classic static solutions and energy of solitons. 
2.1 Static solutions of 4 . 

In the case of static solutions of the ip 4 model the potential is determined by the formulas 

V^) = f(v? 2 -y) 2 , (2.31) 
therefore the equation of motion has the form 

(p"(x) + mV - g^ 3 = 0, (2.32) 

that yields (11.21) in the form 

(tpr = k<p 2 -^-) 2 + 2W. (2.33) 
Its restricted solutions, as it follows from phase plane analysis, exist if 

2 

Tfl 

- — < W < 0. (2.34) 
4g ~ ~ 

The separatrix (W=0) solution of (12.331) is the kink/antikink 

[2 m 
cp = ±y-6tanh(6x)), b = -j= . (2.35) 

While inside the interval (12.341) the equation fl 2 . 3 3 1) is expressed in terms of the elliptic Jacobi 
sinus _ 

/2 m 

-kbsn(bx;k), b = — < k < 1. (2.36) 
9 VI + k 2 

The constant W is given by 

w - -^y'i- < 2 ' 37 ' 

The energy in the Nahm case k = i should be studied separately, because the links (12. 36112. 37p 
between k, b and W is not valid. 

The family of the restricted solutions contains also the constant vacuum ones (W=0) 

m 

* = ± vf ' (2 ' 38) 

After the substitution of (12.351) into (11.241) we obtain the following potential u(x) : 

u(x) = -6b 2 /ch 2 (bx), (2.39) 

with the meaning of the constant b = m/y/2. As a result the two- level reflectionless potential 
of one-dimensional Schrodinger equation —d 2 x + u(x) appears. Eigenvalues and the normalized 
eigenfunctions of which are correspondingly (its numeration is chosen from above to lowercase). 

Ai = — b 2 , ip\{x) = a/36/2 smh(bx) / cosh 2 (bx); 

A 2 = -4b 2 , ^ 2 {x) = v / 36/2cosh(6x). 

In the case of the periodic solution of the phi-in-quadro model the potential has a " cnoidal" 
form 

u{x) = -6k 2 b 2 cn 2 (bx; k) + (5k 2 - l)b 2 . (2.40) 

This potential differs from second Lames equation potential by the constant hence its spectrum 
is two-gap one. 
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2.2 SG model 

Let us briefly describe SG model |35j we integrate the equation (11.11) with the potential (11.51) 
arriving at the first-order differential equation with the parameter W (11.21) . 

(VT = ^-(1 + cos{±M<p)) + 2W. (2.41) 
6g m V 2 

The solutions are restricted and hence have the direct physical relevance, if 

' <W<0, (2.42) 



3g 

that follows from phase plane analysis. If W = the nontrivial solutions are interpreted as 
kink and antikink 

IT IT 
(p(x) = ±\ — arcsin tanh (ma;) (mod&), $ = 2m7rW — , (2.43) 
\l 3g \ 3g 

while at the interval 

-^<W<0, 

the solution of (12.411) yields a periodic function expressed via elliptic Jacobi function. To find 
it one plug ip = ±2my^ arcsin z, then the equation (12.41 j) goes to 

{z'f = m> (l + - , 2 ) (1 - , 2 ) (2.44) 

The solution of (12.441) at the interval (12.421) is given by 

z = ksn(mx; k), (2.45) 

where 

is the module of the elliptic function. Hence 

W =^ 2 -^\ ( 2 .47) 



3^ 

Finally 



/ 2 

ip = ±2mW — arcsin ksn(mx; k) (mod&). (2.48) 
V 3<7 

The class of restricted solutions contains also 

= 0(mo#); W = (2.49) 

and 

93 = ±iim^-{mod<$>)-, W = 0. (2.50) 



Other static solutions are obtained by shifts 

x — > x + x , (p — > <p + $, 

that follows from Klein-Gordon equation invariance and SG equation potential periodicity. 

Let us evaluate the energy of the nontrivial static solutions of both models via the energy 
density definitions 

/oo 
M) 2 /2 + V(<p))dx, (2.51) 
-oo 

for kinks and, in a case of periodic solutions, 

E = 2 I e(x)dx, (2.52) 
Jo 

the constant "1" is the period of the solution. 
For the SG kink: 



16m 

E k = , (2.53) 



and, for a periodic soliton 



]m 2 



E p = — -[(1 - k 2 )K + 2E], (2.54) 



where K(k),E(k) - complete elliptic Legendre integrals. 

2.3 Static solutions of Nahm model as a specific case of cp 4 model 

The solution of the (11.111) . is a particular but specific case of static GL model (m=0, g=2) 

0' 2 = 4 - w \ (2.55) 
A solution of 111. lip is expressed via elliptic functions [7], namely while 

r , (2.56) 



_ W 4 J Q ^(02 _ W 2) (02 + W 2) 



or 



<t> 

dt 1 rzr, dt 



i rt; dt _ l f 

W J ° ^/(t 2 -l)(t 2 + l) W J ° iy/(l-t 2 )(l-i 2 t 2 )' 

Hence the solution is 



that may be tested by direct differentiation in (I2.55p . 
The invariants 

92 = w 4 , 

93 = 0, 

alternatively determine the potential in terms of the Weierstrass function 

,3 



(2.57) 



(2.58) 



w 

z,X)=w + - — m . (2.59) 

p(z; w 4 , 0) — w z /2 
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or, if take into account he relation 



p(z; w 4 ,0) = e 3 + 61 63 a = y /e 1 - e 3 , k=J— — — (2.60) 
sn z {az) V e i — e 3 

one has for e 3 = 0, e x = — e 2 = w 2 /2, the parameters of the solutions in Jacobi terms k — i, a — 
w/V2 

= a sn[a(z - zo)]dn[a{z - go)] 
cn[cr(z - z )J 

This form coincideds with one from [5] . We will use the expression of the potential in terms of 
the solution (jXgTj) 

Vfaix)) = = 60 2 = -6a 2 (1 - cn 2 (az,t)) . (2.62) 



3 The generalized zeta-function via Hermit equation 

In a spirit of Hermit approach, see, e.g. [24] , the function gi(p,x,x) = J exp[—pt]gi,(t, x, x)dt 
is a solution of bilinear equation 

2GG" - (G') 2 - A(u{x) + p)G 2 + 1 = 0, G(p, x) = g L (p, x, x) (3.63) 

which in a case of reflectionless and finite-gap solutions is solved more effectively than ( 14.911) . 
It is possible to cover all necessary classes of solutions of the models ( A,B for SG, C,D for 4 ) 
case and Do for Nahm, via the universal representation by means of polynomials (in p) P, Q 

G(p,x)=P(p,z)/2^/Q{p), (3.64) 

where 

z = sech 2 {bx) 

for kinks A,C, and 

z = cn 2 (bx; k) 

for the periodic B,D. 

b 2 (p(z)(2PP" - (P') 2 ) + p'(z)PP') -{p + u{z))P 2 + Q = 0, (3.65) 
the primes denote derivatives with respect to z, while 

p = { z 2 (l- z), cases A,C; (3.66) 

z(l - z)(l - k 2 + k 2 z), cases D ,B,D. 

b 2 (l-2z),case (A) 
b 2 (2k 2 - 1 - 2k 2 z), case (B) 
u{z) = { 2b 2 (2 -3z), case (C) (3.67) 
b 2 (5k 2 -l-6k 2 z),case (D) 
-Qb 2 (1 - z) case (D ) 
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Substituting (13.691) into (13.651) gives for each power of p — 0, 1, 2 

-2P l {z)-u{z) + q 2 = 0, 
b 2 {2p{z)P^ + p'{z)P[) - Pi - 2u(z)P 1 + qi = 0. (3.68) 
6 2 (p(^)(2P 1 P 1 " - Pf) + p'^P^ - u(z)P? + g = 0. 

respectively. 

Let us start with the cases (A,B). The form of the polynomial is determined from well-known 
facts of the reflectionless potentials theory. The substitution 

P=p + P 1 (z), Q =p 3 + q 2 p 2 + qip + q (3.69) 

into (I3.68P results in 

P 1 = k 2 b 2 z, q 2 = b 2 (2k 2 - 1), q x = b 4 k 2 {k 2 - 1), g = (3.70) 
includes (A) in a sense that for the case k = 1 (13.701) gives 

Pi = b 2 z. (3.71) 

Going to the cases (C,D), generally 

P = p 2 + P 1 (z)p + P 2 (z), Q =p 5 + q 4 p 4 + q 3 p 3 + q 2 p 2 + qip + q , (3.72) 

A substitution of (I3.72p into the Hermit equation splits in the system 

-2Pi -u + g 4 = 0, 
-2P 2 - P 2 - 2uP l + b 2 (2pP{' + p'P[) + q 3 = 0, 
b 2 (p(2P>> + 2P 1 P{' - (P[) 2 ) + pf {PI + P 1 P[)) - 2P 1 P 2 ) - u(2P 2 + Pi) + q 2 = 0, (3.73) 
b 2 {2p{2P['P 2 - P[P' 2 + P X P'{) + p'(P x P' 2 + P[P 2 )) - P| - 2uPxP 2 + gi = 0, 
b 2 (p(2P 2 P% - P 12 ) + p'P 2 P' 2 ) - uP 2 + g = 0. 

The potentials u(z) = c — 6b 2 k 2 z, yields for the case D 

P 1 (z) = az + P, P 2 (z)=a'z 2 + b'z + c', (3.74) 

where 

a = 3k 2 b 2 , 
(3 = 3b 2 , 

a' = I8b*k\ (3.75) 
b' = -3b 2 k 2 (llb 2 k 2 - q 4 - b 2 ), 
d = Q3b 4 k A /A - 5b 2 k 2 qi /2 + b 2 q A /2 - g|/4 + g 3 - 9b 4 k 2 /2 + 36 4 /4, 

are functions only of k, b, where 



q = 0, qi = -27k 2 {l-k 2 ) 2 b 8 , 
q 2 = -9b%k 2 + l)(fc 4 - Ak 2 + 1), g 3 = 36 4 (1 + 9k 2 + fc 4 ), g 4 = 5b 2 (l + k 2 ) 



Finally, 



i=5 



i=i 
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(3.76) 



Q = H(p- Pi ), (3.77) 



where the polynomial Q simple roots pi are ordered so that (0 < k < 1) 

- (2Vl - k 2 + k* + 1 + k 2 )b 2 < -3b 2 < -3k 2 b 2 < < -(2Vl - k 2 + fc 4 - 1 - £; 2 )6 2 . (3.78) 

For the cases (Do? C, D) we begin from that it is the result of substitution of Pi from the first 
equation of (I3.68|) into the second one. Next, the third equation yields 

Q = p(p 2 + (1 - k 2 )b 2 )(p - k 2 b 2 ), (3.79) 

with the simple roots Pi. while in the particular case of 4 (C) 

go = 0, q x = 0, g 2 = 366 6 , g 3 = 33b 4 , g 5 = 10b 2 . 

The arguments in (13.731) are omitted. 

Let us pick up the expressions determining j(p): 



J §th = i^I dx+ ^I {az+ p)dx +i mS {a ' z2 + h ' z + d)dx - (3 - 80) 



t(p) 



The case of Nahm equation yields q^ = 0, q% = —21b 4 , q 2 = 0, qi = 1086 s , go = 0, hence 
PJz) = Sb 2 (z - 1), P 2 = 186V - 36b 4 z. 



i=5 



Q(P) = T[(p-Pi) = p(p + 3b 2 )(-p + 3b 2 )(2V3b 2 - p)(2V3b 2 + p), (3.81) 



i=i 



where the polynomial Q simple roots pi are easily ordered for real b. 
So we need three integrals over the period K. 



J Q dx, J Q zdx, J Q z 2 dx. (3.82) 

Let us go to the variable z, dz = d(cn 2 (x)) = —2cn(x)sn(x)dn(x)dx, a bit more convenient to 
put y = 1 — z, dy = —dz 

/» -V; = | /J ^,.*. 1 ,„ = (3.83) 
/* Sn <(*; fc)<ix = /„> ^?2- iH) = 5*« 2 + * 2 ) A '<*) - 2 <! + *= 2 > B «>- 

hence 

7(p) = 2K(^)^ + a(-™^ + (/3 _ «)2iT(A;)^+ 
^ (a+^W-ga+m*) + (6 > _ 2a')^^ + (a' -b' + d)2K(k))jfa. 

or, finally 



(3.84) 



C(«) = y (jl(p)e pt dp)t s - l dt/T(s); (3.85) 



which is expressed via integrals 



Ji ^ (3-86) 
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by a contour that contains all branch points of the integrands (inverse Laplace transform). Or, 
via 

e^f-Ht = / e\-) s - l d- = — — / e-H^dt = -_±_r(s), (3.87) 



'0 Jo 
Imp < 0, one arrives at 



> V {-P) s ~ 2 Jo (~P) 



In the case of Nahm the mass depends on b: 

P > 2^p(p + 3b 2 )(-p + 3b 2 )(2VSb 2 -p)(2VSb 2 + p) 

The result is given by hyperelliptic integral that can be evaluated numerically or via elliptic 
functions if the symmetry of integrand is taken into account. 

4 The generalized zeta- function as a combination of Baker- 
Achiezer functions. 

Let us consider the Green function defined in the Sec.l by the relation fll,19p . The method we 
use here is based on the technique of finite-gap integration of the KP equation [251 12SI [27] . 

The source part of the problem is one-dimensional version of the heat equation (11.221) with 
the potential u(x) = V"(<p(x)) 

d d 2 \ 

-— — — — — -|- u(x) G (t, x, x ) = 5(t)5(x- x ) . (4.90) 



dt dx 2 

Let us integrate the equation (I4.90p by x over [x — e, Xq + e], hence 



lim e _ /* o °+ e * [j t ~ -^2+ u(xj) G(t, x, x )dx 



(4 91) 

lim^^- M^og^ + dG^o-^o) ] = lim _ o j*o* s {t) 6{x _ Xo)dx = 6{t y 

The function G(t,x,x ) is supposed to be continuous at t > 0. 

Solutions of the equation (14.901) with the zero r.h.s. are expressed in terms of the Baker- 
Achiezer functions ip s (x,t;P) built by the polynomial q(s) = sx + s 2 t of local parameter s 
[25| [27] . Consider a Riemann surface of genus g and a theta function on it. Then there are 
holomorphic differentials u>k, k — 1, ...,<?, normalized as 

u k = 2n5 ik . (4.92) 

The non-special poles divisor on the surface is denoted as D = Pi + ... + P g and the matrix of 
periods is 

B, jk = [ u k , (4.93) 
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that define the general Riemann Theta: 

e(z) = Q(z\B) = exp[(BN ■ N)/2 + (N ■ z\. (4.94) 

Nez g 

The parametrization by the polynomial q(s) yields 

1>.(x,t;P) = C(P)eM J' [dflWx + <B)F>t\) 6(g e ( ^V + g| ^ ' (495) 
where A(P) = j£ uj is Abelian map, C(P) is a constant and the vectors of b-periods 

Ui = I dtt^, Vt = / c/fi (2) (4.96) 



define the argument of the 0-function. The link between the potential u(x) and ip is given by 

u(ar) = -2— Inetetf + tV + D) (4.97) 
ox 2 

which is recognized as general Matveev-Its formula |25j . 

So the parameters of the potential and the Green function arise from (14.96P as the position of 
the cycles bi depends on parameters of the potential via the curve that determine the Riemann 
surface. 

Applying the classical method, one builds the Green function which has the form prescribed 
by (I4.9ip via the linear independent solution ip- s (x,t; P) as: 

G {t, x, x ) = mMx> *; P)^-s{xo, t; P) (4.98) 

that account the unit jumps of the first derivative with respect to x at x = xq and the Wronskian 
M = ipsip'-s ~ ^-s^'s is determined by a normalization at the borders of the period. 

The boundary condition at t = is satisfied if integrate (I4.98P with respect to s along a 
contour C determined by the spectrum of the operator D 1 

d S ^iP s {x,t;P)^ s {x ,t;P). (4.99) 

The zero value at t < is guaranteed by the theorem on Laplace transform [7j. 
The diagonal values of the Green function are given by 

G(t, x, x) = J c § ij}„{x, t; P)^- S (x, t; P) = 

r^QPvr,(T P WO -I- f P - jn\\ Q(xU+tV+A(P)+D) e(xU+tV+A(P.)+D) (4.1UU) 

j UA exp(j Poo a\t -r j Px a\t)) e(xU+tv+D) e{xU+tv+D)) 

P_ corresponds to q- s = — sx + s 2 t. The function j(t) and, next, the zeta-function are imme- 
diately written as fll.2ip . 

The case of SG illustrates the idea in the simplest manner. We start from the Lame equation, 

d 2 \ 

2p(u) * = HV, (4.101) 



du 
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which solutions 

«= g( " + g ) ( ~ g)) exp[c(f-(-g)M. (4.102) 

form a complete set (see e.g. [27J EI])- The linear independent solution is obtained by the 
change of the sign before jp _1 (— H)). 

The equation f)4.10ip is directly connected with one with the cnoidal potential (11.241) 

d 2 \ 

' + 2k 2 - 2k 2 cn(z, kf )* = (4.103) 



dz< 



where z = iK'+u^T^, H= (e x - e 3 ) h + 2e 3 ,k 2 = ffeff , ei = (2- A; 2 ), e 2 = (2/c 2 -l),e 3 = 
— (1 + k 2 ), and v. v. = p(ui),Ui = u,u 2 = u> — to 1 , u 3 = u>'. The link between Weierstrass and 
Jacobi functions 

sn 2 (w, k) = — \ — — , w = zy/ei — e 3 , (4.104) 
p{z) - e 3 

is obviously used. It allows to express a solution of (I4.103P as 



* = ~ , exp[^= ? (p- 1 ((e 1 - e 3 ) + 2e 3 ))], (4.105) 



+ p- 1 ((ei-e 3 )/i + 2e 3 )) p z - iK' 



^ B== Vei - e 3 



e x - e 3 = 3, p-\H) =p,h= P [ P ^ 3 ; u = ^7== ^'( fc/ ) is the complete elliptic integral. The 
parameters are expressed in terms of half-periods of the Weierstrass function, defined by the 
curve with the uniformization (p'(u)) 2 = 4 (jp(-u)) 3 — g2p(u) — g 3 . 

It is convenient to perform numerical calculations using Jacobi theta-function , namely 

a(u) = 2u;exp[^-] ^, 2 ^ , r\ = <;(uj), that leads to 

ggP^gig (^) 

« (fc) 1 ' 

The particular elliptic solution of the heat equation with the potential as in (14. 103[) is the 
product ip(z,t, h) = exp[—ht]^(z). 



5 BA function and "twisted BPS monopoles" 

The BA function (Its-Matveev formula) for SG case is expressed in terms of Jacobi theta 
functions and solves the correspondent Lame equation with ul{x) = 2k 2 (l — cn 2 (x, k)) and the 
spectral parameter h [7J 

exp ^y 2 * +v({p- 1 (3h -2(l + k 2 ))))$(^) 
i; + (x; h) = PL 2 - ™ \ v ^ }}h {2 " J , (5.107) 

where p= p~ l C$h - 2(1 + k 2 )),r] = ((uj),lu = p^ 1 (2 - k 2 ),v = (x-iK')/y/3, K'(k') is the 
complete elliptic integral, p is the Weierstrass function, ((u) is the Weierstrass zeta function 
((' — — p) [TJ-The potential ui{x) of the standard Lame equation differs from one u(x) by the 
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constant factor and shift, namely ul = b 2 u + b 2 , compare with f 1 3 . 6 6 j) . In this section we choose 
for simplicity of formulas 6=1. Hence the spectral parameters are connected by h — p + 1. 
The basic theta function of the representation (15.1071) is defined by 

oo 

$ (w) = d (w\t) = exp[i7r(m 2 r + 2mw)}, (5.108) 

m=— oo 

where exp[mr\ = r], Imr should be positive for the series convergence. The series convergence 
is rapid, therefore the representation (15.1071) is convenient for numeric evaluation of the integrals 
in the zeta formalism. 

The Green function of the spectral Lame problem may be constructed as a product of 
two independent solutions ip + ,ip_ of the spectral equation with the same h: 

, , _ 1 , i(j + {x;h)i;-{xo;h), x < x . . 

g h [x,x )-—{ ^_( x . h ^ + ( XQ . h) x>X(j - ^- iuy J 

The Wronskian factor W is chosen to normalize ( 15. 1 lift so that account the unit jumps of the 
first derivative with respect to x: 

r rdg h (x + e,x ,x ) dg h (x Q -e,x ) f K iin\ 

W dx dx ] = ' (5 - U0) 

The independent solution if)-(x, t; h) may be chosen antisymmetric with respect to the reflection 
x — > —x, e.g defined via "&(w + |). . 

The boundary condition at t = may be account via the integral by spectrum (see also 
(I1.19P ) of the linear independent solutions product 

a(f _ _ a _ f 1 r *+{x,t;h),*-(x ,t;h), x < x 

g{t,x,x )-J —{ 9+ ( XQit . h))i B_( Xit . h)) X>XQ \ dh - ( 5 - m ) 

where ^ + (x,t; h) = exp[—ht]if) + (x; h). 

Finally we integrate the diagonal values of the Green function j Dl (t,x) = J g(t,x,x) from 
(15.1111) by the period of a solution and, after that, one can use the definition (I3.89P of the gen- 
eralized zeta- function via (15.1111) . Integrals dependence on zi are expressed via the Weierstrass 
function of w. It is also known, that 

9(z + B) = exp[---z]9(z), (5.112) 
so the ratio of the theta-functions has the period B. 



6 SG kinks 

The results of the previous sections allow to evaluate corrections to actions for all four (A,B,C,D) 
cases. The results for <p 2 model kinks are well-known, see, e.g. p2], where a table for the 
dimensions d=l,2,3,4 is listed. These results fit the case B, which hence was strictly verified. 

Let us present the formulas for the kinks of the SG model. The substitution of the expres- 
sions from (I3.64ll3.81f3.69f3.71l) yields a divergence of the integral for (13.801) j(p). To regularize 
the Green function let us divide itd Laplace transform into two parts as 

G(p,x) = G c (p) + G k (p,x), (6.113) 
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where the part G c (p) is the Green function diagonal (the solution of (13.631) for a constant 
potential: 

G c = ; . (6.114) 

The kink part is easily constructed via (13. 64113. 81113. 69f3. 711) : 

Gfc = 6W(fe^ (6U5) 
2py V + b 2 

The representation (16.1131) results in two contributions of the quantum corrections to the kink 
(antikink) mass. 



The first one coincides with ( 16.1161) for v 



m 



2 



lDo (t) = (4 7 r)^ r(g + 1 / ) m rf - 1 - 2 -, (6.116) 



and the second one is obtained by the general scheme, namely 



Jk(p) = / G k (p,x)dx. (6.117) 

J — oo 

Plugging (16.1151) into (16.1171) yields 

J k (p) = yL= . (6.118) 

PyP + 

The value of the integral Jf^ sech 2 {bx) = 2/6 is taken into account, the corresponding function 
(ll.2ip or, more exactly (??) will be denoted 7^, the index R label the result of a renormalization 
provided by the division (I6.113p . The function may be found directly from a table Q29J) but 
we would explain the result as an example for further development of the renormalization 
procedure. Note that the cut for the radical \/p + b 2 is made along the Re p -axis from —00 to 
— b 2 , the branch with \Jp + b 2 = iy/\p + b 2 \ on the upper and \/p + b 2 = — iy/\p + b 2 \ on the 
lower bounds of the cut is chosen. After the deformation of the integral contour one has for 
t > 

R1 b f°° exp(£ + b 2 )t 



rt = 1 -i (elm* *- (6.H9) 
^ Jo {£ + b 2 )V£ 

The formula (16.1191) gives an expansion of the resolvent of the operator d t + L by the operator L 
spectrum. The direct substitution of ( 16.1191) into (??) leads to a divergent integrals. However 
in this case a renormalization is not necessary, the integral in the r.h.s. of ( 16. 119ft is expressed 
in terms of the error function 

Erf(z) = 4= / e- r2 dr. 
V 71 " Jo 

namely, after some change of variables 

7 fc (t) = Erf(bVi) = / exp[-b 2 tr 2 ]dr. (6.120) 

Jo 
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Figure 1: d=l,2,3 dependence of the derivative v , — - to mass on the parameter m 

The same result gives [29]. The Mellin transform of this representation gives the following 
expression of zeta function of the operator L. 

__2&- 2 *r( s + i/2) /* b-*°r( s + 1/2) 



£?(s) = = v 7 7 / r-^dr = — = v " '7 . (6.121) 

M ; v 7 ^ r(s) y r( s + i) 1 ; 

The integral in (16.1211) converged only at -Res < but gives the analytical continuation for all 
s G C excluding the poles of T(s + 1). 

Substituting the result f)6.12ip into one arising from (11.251) with account of (I6.116P yields 

C f ( S ) = -4(47r)fm' j - 1 - 2 % r(s + 1 ~, rf/2 \ . (6.122) 
SLV ; v ; (2s + 1 - d)r(s) v ; 

The condition b = m is taken into account. Differentiation of (16.1221) by s at the point s = 
gives the desired correction 



<CfW) 

2 ds 

-2(4n)lm d ~ 2s - 1 ^g^ ((d - 2s - 1) (Psi (-id + s + l) - 2 lnm - Psi (s)) + 2) 

(6.123) 



next we plot the dependence of the correction ^L^i Q n m (Fig [p. 



We would remind about the choice of the constant g=2 in the last sections. 

7 Conclusion 

The integrals in the elliptic case are evaluated numerically by means of rapidly converging series 
for theta-functions. 
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It is known [35] that it is possible to form a periodic solution of a soliton model by a shift 
operation in the complex plane of the soliton (kink) parameter. 

In [37] the authors study the diffusion of kinks. The Ginzburg-Landau (GL) model is very 
popular in different aspects of solid state physics, e.g. for magnetics [SB]. Some recent papers 
open new field for applications [321 [161 EH- 

The investigation of dislocations dynamics by means of FK model gives a direct possibil- 
ity to check quantum soliton effects separating kink and elliptic solitons contribution via the 
energy dependence on parameters [37J. It is interesting to incorporate our results in a real crys- 
tal thermodynamics via statistical physics approach [38J or a direct echo response evaluation 
with correspondent measurements [39]. In a review [38] a quantization contribution is already 
discussed. There also direct simulations of kink-antikink pairs [40] . 

7.0.1 Acknowledgment 

I thank G. Kwiatkowski for discussions and assistance in plotting of the SG kinks mass depen- 
dence on m. 
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